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Abstract 



O ' Structural transformations, electronic spectra and ballistic transport in 

pulled gold nanowires are investigated with ab initio simulations, and corre- 
lated with recent measurements. Strain- induced yield of an initial double- 
strand wire results first in formation of a bent-chain which transforms upon 
further elongation to a linear atomic chain exhibiting dimerized atomic con- 
figurations. These structures are stabilized by directional local bonding with 
spd-hybridization. The conductance of the initial double-stranded contact is 
close to 2(2e 2 //i) = 2g$ and it drops sharply to l<?o during the transformation 

q . to a single chain, exhibiting subsequently a ~lgo plateau extending over an 

elongation well above typical Au-Au distances. 
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Generation of nanometer-scale crystalline wires (NWs) through elongation of materials' 
contacts had been predicted via early molecular dynamics simulations jl] , and their mechan- 
ical, structural, and electrical properties have been the subject of intensive investigations 
owing to increasing basic and technological interest in such NWs. Most recently Q 
high-resolution electron microscope (HRTEM) images, recorded during retraction of a tip 
from a gold surface, portrayed the development of a sequence of crystalline NW structures 
made of parallel atomic strands oriented along the NW axis, with the number of strands 
decreasing one-by-one culminating in an one-atom wide and four-atom long linear chain of 
Au atoms [|3],^| . Furthermore, the simultaneously measured conductance revealed [Q a stair- 
case of quantized values (close to integer multiples of the conductance quantum g Q = 2e 2 /h) 
with the disappearance of each atomic strand, coinciding at the final pulling stages with 
an approximate unit (go) drop in the measured conductance. These experiments stimulated 
recent electronic structure calculations pertaining to single-chains of gold atoms. [|5|-§] 

Aiming at elucidation of the atomic-scale structural evolution and transformation mech- 
anisms, and of bonding and transport characteristics at the ultimate stages of elongation 
of gold NWs, we report here results from large-scale ah initio density functional simula- 
tions P,|10|| in conjunction with conductance calculations |fLl |. Our findings reveal that the 



formation mechanism of a four-atom long single chain NW, generated through straining 
of a double-strand (ds) wire (see bottom left configuration in Fig. la), involves stress ac- 
cumulation followed by eventual yield of one of the atomic strands, with the consequent 
atomic rearrangement resulting first in a (lower-energy) bent-chain (be) configuration (top 
left configuration in Fig. la). The yield event is accompanied by a significant decrease in 
the pulling force and a sharp drop in the conductance to lgo (see top in Fig. la). Further 
separation of the two electrodes results in a transformation to a linear-chain (lc) struc- 
ture (see right configuration in Fig. la), with continued elongation leading to formation of 
dimerized configurations. Initial dimerization develops between the inner atoms of the wire 
(i.e. inner-dimerization, id), transforming at a later stage to one where each of the inner 
atoms pairs with an end-atom of the chain (i.e. end-atom dimerization, ed). During the 
ds— >bc— dc evolution the conductance exhibits a ~lgo extended (4 A - 4.5 A long) plateau 
(see Fig. la). The optimal structures determined here are stabilized by directional local 
bonding involving s,p and d hybridization near the Fermi energy, and they are found to be 
energetically favorable to certain suggested alternative structural models. [[F] 

We start from a relaxed ds wire configuration with the distance between the two out- 
ermost layers of the opposing Au electrodes held at L = 11.5 A; each of the electrodes is 
composed of 29 atoms arranged in 3 layers with face-centered-cubic (110) stacking in the 
z-direction (along the wire) and exposing {100} and {111} side facets. This crystallographic 
orientation, and the subsequent elongation process and conductance characteristics discussed 
below, correspond to experiments discussed in Ref. 3 and described in Fig. 3 of that pa- 
per, where HRTEM and conductance measurements pertaining to gold nanowires created 
via controlled elongation of the contact between a tip and a (110) facet of a gold island 
are displayed. Note that this experiment is distinctly different from the one corresponding 
to Fig. 4 of Ref. 3, which pertains to a wire formed through (spontaneous) thinning of a 
"nanobridge" connecting across a hole formed in a gold film by intense electron bombard- 
ment (in particular, unusually large interatomic distances of up to 3.5 A - 4 A were reported 
in Ref. 3 for the latter case, and not for the tip-pulling experiments which are the focus of 



our study). 

In this initial configuration the 2-atom strands are parallel to each other, with al- 
most equal inter- and intra-strand separations, d intra (l', 2') = d intra (l,2) = 2.66 A and 
dinteri^-' ', 1) = d in t er (2',2) = 2.67 A (see atom numbering in Fig. la), and the distance be- 
tween the topmost (interfacial) facets of the opposing electrodes is d e i- e i = 6.02 A. Note that 
all the distances involving atoms bridging the electrodes are smaller than the nn interatomic 
distance in bulk gold (2.885 A), but larger than the bond-length of the free Au2 dimer (2.48 
A). The calculated conductance [JXTJ of this configuration is 1.79g jl2 |- 



Elongation of the wire is simulated through increasing the separation between the outer- 
most layers of the opposing electrodes (where the atoms are held at their fee lattice positions 
with a lattice constant of bulk gold, 4.08 A), followed by full relaxation of the system after 
each elongation step. The total elongation simulated here is AL = 6.5 A (using elongation 
increments of 0.25 A< 6L < 1 A, see Fig. la), and we found that structural changes in 
response to increments in L involve mainly atoms forming the connective wire bridging the 
electrodes. 

ds^bc transformation (L < 13.3 A). Initial elongation results in higher-energy strained 
configurations with the pulling forces rising to ~ 4 nN (corresponding to the elongation 
increment from L = 12.3 A to 12.6 A, see Fig. la), and it is accompanied by a gradual 
(small) decrease of the conductance to G = 1.68go i n the highly strained ds configuration 
shown at the bottom left in Fig. la (corresponding to L = 12.6 A with d e i_ e i = 6.95 A). 

For all the ds configurations the ballistic electron transport involves mainly two conduc- 
tance eigenchannels (CCs), with the total conductance G = go J2n \ T n\ 2 , where < \r n \ 2 < 1 
is the transmission probability of the nth eigenchannel; for the initial ds configuration 
|r x | 2 = 0.99, |r 2 | 2 = 0.76, and for the highly-strained one \r x \ 2 = 0.97, |r 2 | 2 = 0.70. The 
CCs are determined by quantization of the electron motion transverse to the propagation 
direction (z), due to the lateral confinement by the effective wire potential Jl0|JTl|. Each of 



the transmitted CCs is found to be delocalized over the two atomic strands with the first 
CC being nodeless, and the second having a nodal plane normal to the plane defined by 
the two strands and bisecting the interstrand Au-Au bonds (see Fig. 2a); the delocalized 
nature of the CCs contrasts the conjecture || that in such multi-strand NWs each of the 
channels is associated with (and is confined about) an individual atomic strand implying 
two independent CCs connected in parallel to the electrodes. 

The ds— >bc elongation stage culminates in breaking of one of the strands, accompanied 
by displacements of each of the end atoms of the broken strand (1' and 2' in Fig. la) to 
the center of the underlying (110) rectangular facet of the corresponding electrode, yielding 
a lower-energy be configuration, with a concomitant sharp drop in the conductance to 1 go 
(see Fig. la) involving a single CC. In the relaxed be configuration (top left structure in Fig. 
la) the end-to-end wire length l ee = d(l', 2') = 5.12 A, the distance between an interior atom 
to an nn end-atom of the wire d ie = d(l, 1') = d(2, 2') = 2.62 A, the bond length between 
the internal atoms of the wire da = d(l, 2) = 2.56 A, and the bond angle 1(1', 1, 2) = 119°. 

bc^lc transformation (13.3 A < L < 16.3 A). During subsequent elongation, requiring 
initially a small pulling force (~0.8 nN, corresponding to the interval between L = 13.3 A and 
L = 13.8 A), the be wire straightens gradually resulting eventually (see right configuration 
in Fig. la) in a linear chain at L = 16.3 A |13[ . The be— dc transformation spans an 



elongation range AL(bc — > Ic) = 3.0 A, during which the conductance decreases only by 



about 10% (from 1 go for the initial be configuration to 0.87 go for the lc one at L — 16.3 

A). 

In light of a recent report pertaining to alternative structures for single-chain gold 
nanowires we have explored such configurational isomers during the be stage, with two of 
them shown in Fig. lb; a planar 4-atom rhombohedral (rh) arrangement at L = 13.3 A and 
a "zig-zag" (zz) structure at L = 14.8 A. Energetically, these (relaxed) structural isomers 
are found to be local minima with higher energies (and slightly lower conductance) than the 
corresponding be configurations (see Fig. la). The zz isomer and be structure (as well as the 
lc) become approximately degenerate at L = 15.8 A, and upon further elongation (i.e. for 
L = 16.3 A) both the zz and the be convert to a linear chain. Note that although all these 
isomers (rh, zz, be, lc) are found to be local minima for a free AU4 cluster, the energetic 
ordering of structural isomers observed here for the four-atom NW cannot be deduced from 
that of the free cluster |14| . 

Ic^breaking (L > 16.3 A). With continued pulling the lc wire shows a tendency toward 
dimerization. At first (for 16.3 < L < 17.3 A) the only stable structures correspond to 
internal dimerization (id) of the wire, with da < di e [[HJ. At the start of this elongation 
interval the conductance decreases only slightly, with a marked drop at L = 17.3 A. However, 
starting at L = 17.3 A an energetically competitive dimerization mode of the lc wire emerges 
(being essentially degenerate with the id structure at 17.3 A) where da > d ie |1|, with the 
conductance of this end-dimerized (ed) structure being lg . We note here (see Fig. la) 
the markedly higher conductance of the ed structure (G = lg ) compared to that of the 
corresponding id chain [[17J. Such end-dimerized structures || are energetically favored for 
the rest of the elongation process (that is for L > 17.3 A) ||18|| . The force required for 
elongation during the end-dimerization stage is rather small (i.e. ~0.3 nN in the interval 
17.6 A < L < 18 A), and the conductance of the wire decreases sharply as eventual breaking 
is approached. 

Throughout the structural evolution the electronic spectrum in the wire region is char- 
acterized by a dominant contribution of atomic d orbitals in the energy range Ep — 4.5 eV < 
E < Ep — 1 eV, with s, p, and spd hybrids contributing at energies above and below this 
range. Examination of the local densities of states reveals that the main variations in the 
spectrum in response to the mechanical elongation occur predominantly in the interelectrode 
spatial region (the wire atoms, and the electrode atoms bound directly to them), involving 
mainly states in the vicinity of Ep. 

To illustrate the nature of these states and their high sensitivity to structural variations 
we examine, for selected atomic configurations developed during the elongation process, 
orbital images (Fig. 2(b-d)) and their angular-momentum decompositions about the wire 
atoms, corresponding to electronic states near Ep (i.e. highest occupied molecular orbitals, 
HOMO(j), with the index j = —1, —2 corresponding to the next-to-highest occupied orbitals 
in descending order). We observe that for the ds structure (see Fig. 2b and the configuration 
shown at the bottom left in Fig. la) the HOMO is of predominant d xz character {p x ' 2 d^:J 6 , 
here and in the following only the dominant m-components of p-and <i-characters are given) 
and the HOMO(-l) is of strong sp z character (s ' u p z 65 d° z ^ 2 ) on the strand atoms. These 
states evolve after the ds^bc transformation (see top left configuration in Fig. la) into 
a HOMO with s 0A5 p° z 50 (p x 25 p° z 25 d J 2 ) on inner- (end-) atoms, and a HOMO(-l) with an 
almost pure (tilted) d z 2 character on the inner-atom (p x ' 03 d^i) and a (s ' 5 p x m d®f 2 2 ) hybrid 



on the end-atoms of the bent-chain (Fig. 2c). 

The structural evolution during the be— dc stage and the formation of dimerized config- 
urations of the lc are accompanied by reordering of the orbitals near Ep with the HOMO 
orbital on the inner atoms of the wire acquiring a predominant ip x m d° x ^) d xz character 
(see e.g. Fig. 2d corresponding to the ed configuration); as evident the d xz components 
on the two inner atoms of the ed wire combine locally in an antibonding manner, while in 
the (lower-energy) HOMO(-l) orbital (not shown) they are in a local bonding orientation 
with respect to each other. The HOMO (-2) orbital (Fig. 2d) of the ed configuration is 
characterized by a strong s character on the inner-atoms (s°' 86 p° z u d° z '2 3 ) combining to form 
a a-like local bond, and an sd hybridization on the end-atoms of the chain {s QM p°^ Q& dP^_ 2 )- 

The above hybridization patterns and in particular the strong contribution from d- 
orbitals to directional covalent-like bonding characteristics in the (interelectrode) wire re- 
gion, underlie and facilitate the structural transformation mechanism described by us, via 
stabilization of the sequence of strained atomic configurations which evolve in the elongation 
process. Consequently, formation of such (several atom long) single-chain nanostructures, 
reflected in an extended conductance plateau (with close to unit conductance) before break- 
ing 0] may not be found for materials where the above bonding patterns do not occur; 
indeed extended single-chains were not observed (theoretically ([HJ and experimentally [HI) 
for a simple metal contact (e.g. sodium). 

Finally, we remark that the investigations presented here, which deepen our insights into 
the structure, formation mechanism, and conductance of atomic double-strand and single- 
strand nanocontacts between gold electrodes, pertain to the type of tip-pulling experiments 
reported in the first part of Ref . 3 (see Fig. 3 in Ref . 3) , and they are not aimed at explaining 
the extraordinary stability of monoatomic thick wires bridging two sections of a suspended 
thin gold film, where interatomic distances up to 3.5 A - 4 A were reported from HRTEM 
images (see the second part of Ref. 3, particularly Fig. 4 therein). Indeed, the formation of 
wires in the latter experiments involves an initial intense electron bombardment to perforate 
holes in the film and subsequent thinning of the formed " nanobridges" driven by stress and 
strain relaxation processes; additionally, the nature of such experiments differs greatly from 
the controlled tip-pulling experiments, on which we focus here. We remark that none of 
the theoretical studies reported to date [|5|-|| have been able to consistently explain the 
abovementioned anomalously long interatomic distances in the nanobridges |20|| , indicating 



that new theoretical considerations are needed, including examinations of possible effects 
caused by light impurity atoms, which while affecting the structure of the nanowire may not 
have enough contrast to be visible in HRTEM images. ^1] 
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FIGURES 

FIG. 1. (a) Energies (AE, in unit of eV plotted at the bottom, relative to the initial 
ds configuration at L — 11.5 A) and conductances (G, in unit of go? plotted at the top) 
of Au NWs versus the distance between the outermost layers of the opposing electrodes. 
L; the double-strand (ds, squares), bent-chain (be, triangles), and linear-chain (lc, circles) 
intervals are marked at the top. The values for the rhombic (rh) isomer at L — 13.3 
A and for the zig-zag (zz) isomers at L — 14.8 A and 15.8 A are marked by crosses, 
and those for the end-dimer configurations at L — 17.3 A, 17.6 A and 18 A are depicted 
by stars. The elongation force between consecutive structures may be estimated by the 
corresponding slope AE/AL. The atomic configurations (with atom indices) shown as 
insets to (a) correspond to the ds (left, bottom), be (left, top), and lc (right) structures at 
L = 12.6 A, 13.3 A and 16.3 A respectively, (b) atomic configurations for the rh isomer at 
L = 13.3 A (d(l, 2) = 2.76 A, d(2, 3) = 2.59 A, and Z(l, 2, 4) = 124°), and for the zz isomer 
at L = 14.8 A (d(l, 2) = d(3, 4) = 2.55 A, d(2, 3) = 2.52 A and Z(l, 2, 3) = 119°). 

FIG. 2 (color) (a) Iso-surface images of the magnitudes of the first (0.97go, top) and 
second (0.7 go, bottom) transmitted eigenchannels for the highly-strained ds configuration 
(at L = 12.6 A), (b and c) HOMO (top) and HOMO(-l) (bottom) orbitals corresponding 
to the highly-strained ds NW (in b) and to the be configuration at L = 13.3 A (in c). (d) 
HOMO (top) and HOMO(-2) (bottom) orbitals for the ed configuration at L — 17.3 A. In 
(a-d) Au atoms are depicted by red spheres. In (a) we display the entire atomic system 
including parts of the jellium-slab continuations of the electrodes (used in the conductance 
calculations), and in (b-d) we focus on the vicinity of the nanowires. The viewing orientation 
was chosen in each case to enhance visualization of the orbital characters on the wire atoms. 
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